Numerical study of magnetic field induced ordering in BaCuSi2 06 and related systems 
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Thermodynamics of spin dimer system BaCuSi2 06 is studied using a quantum Monte Carlo 
calculation (QMC) and a bond-operator mean field theory. We propose that a new type of boson, 
which, rather than being hard-core, allows up to two occupancy at each site, is responsible for 
the Bose Einstein condensation of field induced ordering. Its superfluid density is identified as the 
square of the in-plane staggered magnetization ra^y in the ordered phase. We also compare our 
QMC result of the spin Ideisenberg model to those predicted by mean field theory as well as by the 
simple hard core boson model for both large and small intra-dimer coupling J. The asymmetry of 
the phase diagram of mxy{h) of small coupling J in related systems such as NiCl2-4SC(NH2)2 is 
explained with our new boson operator. 

PACS numbers: 75.10.Jm, 75.40.Cx, 75.40.Gb, 75.40.Mg 
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I. INTRODUCTION 

Recent experiment on BaCuSi206 under strong exter- 
nal magnetic field observed a A-like transition of the spe- 
cific heat capacitjii. This is among other spin dimer sys- 
tems, such as KCuCls and TlCuCla, that exhibit quan- 
tum phase transition from spin liquid to a magnetically 
ordered state with increasing magnetic field. The or- 
dered phase is characterized by an uniform magnetization 
accompanied by a long range staggered magnetic order 
perpendicular to the field. The staggered magnetization 
vanishes as temperature rises above the transition tem- 
perature Tc ■ This has been interpreted as a phenomenon 
of Bose-Einstein Condensation (BEC) of the magnetic 
field induced Sz — +1 triplets, which, when ignoring 
higher energy interactions, behave like repulsive hard- 
core bosonaii. The external field now plays the role of 
chemical potential and controls the number of triplets. 
Measurements of Tc near critical field he are reported 
to have the deduced critical exponent v — 0.63(3)^, 
in good agreement with the predicted value v = ^ oi 
Bose-Einstein condensation. Numerical calculations also 
found the critical exponent approaching the expected as 
h — > hc^'^ for 3D spin dimer systems, support the notion 
of BEC. 

Previous numerical calculation on the hard-core boson 
model found a similar transition as observed!'^. Instead, 
in this paper, we employ the original Heisenberg model to 
study the temperature and field dependence of thermo- 
dynamic quantities. In contrast to the hard-core model, 
we introduce a new type of boson, which allows up to 
two occupancy at each site, and show that its condensate 
density is identical to the staggered magnetization miy 
in the ordered phase with the global phase corresponds 
to the direction of rrixy The temperature and magnetic 
dependance of m^^ is explained in the context of Bose- 
Einstein condensation. We also extend to the case with 
smaller inter-plane coupling in which higher energy states 
are no longer negligible in the ordered phase and leads to 
the asymmetry of on both mz{h) and m'^y{h). This con- 
dition is relevant to other dimers material like TlCuCla 



and = 1 system NiCl2-4SC(NH2] 
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II. MODEL HAMILTONIAN 

In BaCuSi2 06, while Cu^+ ions arrange itself in lay- 
ers of square lattice, every two Cu-Si-0 layers are sepa- 
rated by planes of Ba^+ ions, and the distance between 
adjacent bilayers, 7.043 , is much larger than the bi- 
layer separation, 2.73 A. Consequently, the spins of Cu^+ 
ions form dimers between bilayers and interact weakly in 
plane and out of plane. We study the antiferromagnetic 
Heisenberg model in an external field hz written as: 



-|-J2^Si,j • S2,.,+£ - fj.B9oh^S^ .^, (1) 



where a = 1,2 denotes the layer index and the second 
summation refer to summing over all nearest neighbors in 
the xy plane for both types of layers. The exchange cou- 
pling constants are taken as J = 4.45 meV, Ji=0.58 meV 
and J2=0.116 meV, which is provided by ref. Q whose 
QMC based on the hard-core boson model using these pa- 
rameters yield the best fit to the experimental result. The 
magnetic field h yields the Zeeman energy in the Hamil- 
tonian and the gyromagnetic constant go — 2.306 in the 
case of BaCuSi206. We first focus on the condition that 
J ^ Jl, J2 so that the zero field ground state only com- 
posed of inter-plane singlet dimer states. Finite energy is 
needed to excite a sing let state |s) = (| n)-| iT))/V2to 
a triplet state |t+) = | tt)- Increasing the magnetic field 
h will reduce the energy gap which closes up at a critical 
/ic ~ 24 T such that fJ-sgohc is of order J. For h > h^., a 
new magnetic order with staggered magnetization in the 
xy plane emerges, with the ground state a linear com- 
bination of |s) and in each dimer. The breaking of 
rotational symmetry in the xy plane implies the existence 
of a Goldstone mode at Q=(7r,7r,7r). Raising the field h 
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further leads to a finite transition temperature Tc which 
reaches a maximum before falls down again to zero at a 
saturation field hs ■ 

From another point of view, the quantum phase tran- 
sition at he can be considered as BEC of bosons fot = t^j^ 
on the vacuum composed of \s). Since no more than one 
triplet is allowed in the same dimer, the triplet dimers 
6^, with S=-|-l, behave as hard-core bosons. These 
bosons interact repulsively because Ji and J2 favor anti- 
ferromagnetic couplings of neighboring dimers and hence 
neighboring bosons are avoided. One can reduce H^p to 
an effective Hamiltonian Hif* of hard core bosons by pro- 
jecting out higher energy states |io) = (| Ti) + | iT))/V2 
and = I ID, 

Hb = tiY^ {b\bj + h.c.) + Vi^ niUj 

i i 

■i 

where ti = Vi = Ji/2, t2 = V2 = J2/4 and ni = b\bi. 
The chemical potential, /i = J — i-iB9oh, now depends 
on the magnetic field h and only when h > he, bl finite 
concentration of bosons arises and condensate forms in 
the ground state. Note that this Hamiltonian preserves a 
particle-hole symmetry that the phase diagram of Tc{h) 
(or no{h), the superfluid density), must be symmetric 
about a /i„i where Tc is a maximum at uq — 1/2. Nu- 
merical calculations based on this hard-core boson model 
reproduce thermodynamical quantities that agree well 
with the experimental results and support this boson 
picturei. Matsumoto et al. pointed out that, however, 
when constructing a consistent mean field theorjiS, the 
high energy state should be included. This is reflect- 
ing the process of annihilating and at neighbor- 
ing dimers while creating a \s) and vise versa. This pro- 
cess becomes important when the inter-dimer coupling 
Ji and J2 are comparable to intra-dimer coupling J and 
when h is not much larger than he- Therefore, to provide 
a consistent and general description for all J and h, we 
introduce the appropriate boson operator b^: 



b! 



i-iy 

V2 



(3) 



with its operations on the spin states 





bi\t+)i = \s)i 


^\t+)^ = 


bi\s)i = |t_)i 


bl\t-h = \sh 


k\t.), = 



for even i and an extra negative sign for odd i. Here 
is redefined as — | tT)ij and the first and second arrows 
denote the spins of layer 1 and layer 2 respectively. These 
yield the expectation values {s\b''b\s) = {t+\b''b\t+) = 1 



and = 0. Unlike 6t, 6t operates on a Hilbert 

space of vacuum |0) = while |s) and correspond 
to single and double occupied states of b respectively. \to) 
is again decoupled from the other states and not included 
in the Hilbert space. Although occupation greater than 
two at each lattice site is again prohibited, b boson is no 
longer the simple hard core boson described by b. Note 
that the vacuum is higher in energy and the lowest energy 
state |s) is fully occupied when h < he- 

The modified effective Hamiltonian is now written as: 

i i 

+ti^{b^J}-j + h.c.) + Vi^ fi.fij 

(id) (id) 

+t2j2(brbi+z + h.c.) + V2J2 n^nt+z, (4) 

i i 

where e = — J, fi = —HB9oh and = blbi — bib] with 
ti^2 and Vi^2 unchanged. The potential energy terms arise 
from the term SfSj of the original spin Hamiltonian Hgp 
which is attractive for and t- but is repulsive for the 
same kind of triplets. This Hamiltonian now loses the 
symmetry between |s) and due to the presence of 
\t-). For he < h < kg, it is the condensation of this boson 
bi that gives the quantum phase transition observed. 



III. MEAN FIELD THEORY 



A mean field condensate ground state of the effective 
Hamiltonian Hj, is written asi 



1*0) = ll[ubl + vi-l)\fe^Obf+ge-^0))\0),, 

i 

= n ("l^"^)^ + vi-l)%fe^'\U). + ge-^'\t^),)X5) 



1 and all parameters 



where u'^ + v'^ — 1 and -I- g'^ 
are real. The last equation is exactly the same as the 
one taken for mean field condensate using bond opera- 
tor representation^ in which 9 is chosen to be zero. The 
global phase 9 corresponds to the angle of rotation in 
the xy-plane whose presence, due to the rotational in- 
variance, should not change the energy of the system. 
t+ and i_ undergo the transformation e^^t^ and 

t- — > e~*^t_ when the x and y axes are rotated by an 
angle 9. As shown below, this phase 9 also specifics the 
orientation of the in-plane staggered magnetization. 

Remarkably, the ground state |^o) has two finite ex- 
pectation values that correspond to two order parame- 
ters: 



(b) = uv{f + g)e'' 

(62) ^ y^fge^^<^. 



6Zi 
oe 



(6) 
(7) 



The first one is the usual order parameter expected from 
a BEC, while the second additional one is originated from 
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the fact that b allows up to two bosons occupied at each 
site i. These expectation values can be related to the 
staggered magnetization and spin-spin correlation func- 
tion when Eq. |31is applied. 

In the spin language, the staggered magnetization is 
taken as the order parameter for the ordered phase and 
defined as 



2 2,2 



x,y. 



As mentioned above, one can easily show that rrixy is 
identical to the magnitude of superfluid order parameter 
bo- Using the definition of fet from Eq. O one obtains 



2N 
— bo cos 9. 
Similarly, my = bo sin 6 and consequently, 
rrixy = foo = uv{f + .g), 



(8) 



and the phase 9 specifies the direction of m^y The real- 
ization of m^y — bo allows us to compute the staggered 
magnetization directly using QMC. 

The uniform magnetization along z axis can also 
be written as 

N 

.,2/j-2 „2\ 



Wh-kbi) 



(9) 



which is simply the difference between the number of 
and of \tJ). 

We stress that the distinction between 6^ and fe^ origi- 
nates from the consideration of , which we take as the 
vacuum here. If the \t^) state is much higher in energy 
and being ignored, 6^, as g — > 0, is just the hard core 
boson 6^ described in ref. 1. This approximation be- 
comes exact when J — s- oo, but will be insufficient when 

J ~ Jl, J2- 

A complete description of the static and dynamical 
properties of BaCuSi206 under apphed field can be ob- 
tained by performing a mean field analysis using the bond 
operator representations. The same method has been 
successfully applied to related compounds of KCuCls and 
TlCuCls by Matsumoto et alS. 

Under bond operator representations, each spin opera- 
tor is replaced by boson operators, one singlet and three 
triplet operators, that operate on a inter-plane dimer 
bond and the original spin Hamiltonian is then trans- 
formed to a Hamiltonian of interacting bosons. These 
are hard core bosons because of the constraint 1 — 

s + t+^H- -I- t\to + t^_t-. Two unitary transformations 
are performed to mix the bond operators in the ordered 
phase to give the appropriate condensate a = {ai)av, 



For convenience we set 9 = 0. This ground state conden- 
sate is in fact identical to |^'o)- The mean field approach 
proceeds as usual by taking the uniform field a and 

minimizing the energy to obtain a set of self-consistent 
equations. The parameters u, v, f and g are now de- 
termined by the self-consistent equations. The particle 
number constraint does pose a problem on the decoupling 
of operators in the above procedure. Naively employing 
a Lagrange multiplier to account for the constraint would 
not close the energy gap in the ordered phase. In ref. 
a Holstein-Primakoff approximation is taken instead by 
assuming the contribution of energy modes other than 
is small in the ground state. This approach successfully 
recovered the Goldstone mode with expected features for 
the ground state. We refer the readers to ref. for the 
details and only present our result here. 



IV. QUANTUM MONTE CARLO 
CALCULATION 



The aim of quantum Monte Carlo simulations is 
twofold: to test the validity of the proposed mean field 
theory and to justify the hard-core boson model. We 
stress that the simulations are performed on the original 
spin Heisenberg model rather than the boson model and 
any observed phase transition is a pure consequence of 
spin exchange coupling of Hgp. 

Having been successfully applied to spin systems in 
external magnetic field in the last decade and demon- 
strated its advantages over the standard worldline ap- 
proach, the stochastic series expansion (SSE)- is the 
method of choice for our problem. The algorithm of mea- 
suring Green's function introduced by Dorneich et al^ 
allows us to compute the spin correlations like (S^Sj^). 
In this scheme, however, since spin must be conserved 
in the ensemble, direct measurement of staggered mag- 
netization nixy, which requires the knowledge of (S^) 
and {Sl[), is impossible. Instead, in our calculation, we 
compute the superfluid density 



no 



(10) 



miy, and 



usi + ve 



In MET, Eq. inland Eq. |H|lead to no = b"^ 

to which our simulated no can be compared. Using the 
definition of 6, Eq. 13 the computation of no involves only 
measuring spin correlations which can be done easily with 
Dorneich's method. Therefore a direct comparison of the 
calculated order parameter by QMC to those predicted 
by mean field approximation as well as experimental data 
is possible. 

All simulations are performed on cubic lattice of size 
12x12x12 typically with 4x10^ update measurement cy- 
cles for each data point, except the measurement of spe- 
cific heat which requires up 1.6x10^ update cycles to 
achieve the acceptable accuracy. 
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FIG. 1: (color online) QMC result of (a) specific heat c„ and 
(b) spin susceptibility x vs. temperature for different mag- 
netic fields. Data for different h are shifted for clarity. 



V. THERMODYNAMICS 

The specific heat and susceptibility at finite tempera- 
ture around the transition have been experimentally mea- 
sured by Jaime et a/i for various magnetic fields. The 
sharp peak of specific heat at Tc has a familiar A shape 
similar to the one found in liquid helium He"*. Our QMC 
calculations based on the spin Hamiltonian Hsp repro- 
duce this result for different h as shown in Fig. ^a). 
A peak of the specific heat c„ signals the BEC devel- 
ops when h > he {he ~ 24 T) and the peak grows with 
increasing Tc as h is increased. The calculated values of 
Tc approximately agree with the experimental result. For 
/i = 37 T, our calculated and the experimental Tc are ~ 4 
K and ~ 3.75 K respectively-'^. Remind that we adopted 
the parameters J, Ji and J2 used in ref. Q in which 
the hard-core boson Hamiltonian Hb was used to repro- 
duce the experimental result. Therefore the discrepancy 
on Tc between QMC and experiment is a consequence of 
ignoring higher energy states in H},. One can certainly 
produce a better fit to the experimental data if different 
couplings than those obtained by ref. jll are used. 

The calculated spin susceptibility Xs in Fig. ^h) also 
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FIG. 2: (color online) (a) Low temperature regime of uniform 
magnetization per dimer for different magnetic fields that 
exhibits BEC. (b) The whole temperature range of for 
field /i = 37 T. (c) Magnetic field dependence of obtained 
from both QMC (T=fK) and MFT together with the number 
density of s and t+ states. 



shows a peak at Tc- The formation of boson condensate 
below Tc freezes the spins in the condensate and therefore 
reduces the response to the external field that leads to 
a drop in Xs- Unlike the specific heat, the transition 
peak reduces when h increases and becomes rather flat 
at h hjfi . 

The phase transition is also visible in the low temper- 
ature regime of uniform magnetization in which a mini- 
mum (maximum) is reached at Tc for h < hm (h > hm) 
as displayed in Fig. Ufa). For h < h^, the low- lying en- 
ergy states have the \s) state dominating and so as T ap- 
proach Tc from above, ruz, which is approximately equal 
to the number of t+ states, reduces. However, below Tc, 
due to the formation of superfluid condensate, number 
of states increases as the fraction of condensate grows 
and therefore raises the m^. This feature is also found 
in the attractive Hubbard model of electrons^^ in which 
the double occupancy, i.e. the number of on-site Cooper 
pairs, increases as T approaches from Tc- 

The role of s and i+ states exchanges for h > hm and 
a maximum of m, is resulted instead. The same tem- 
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FIG. 3: (color online) The slope of mz{h), i.e. the suscepti- 
bility, at T = IK. 



perature dependence of magnetization was also observed 
in another dimer system TlCuClj^^ experimentally and 
a Hartee-Fork calculationi^ on hard core dilute magnons 
explained the qualitative feature of rriz (T) around Tc of 
the proposed BEC. QMC simulations on different cou- 
pling strengtb™ has shown that this is a general feature 
of 3D coupled dimer system while the extrema are miss- 
ing in non-interacting dimers. 

Since the total uniform magnetization is simply the 
difference between the number of \t^) state and of 
further increase of temperature (to T ~ 15 K for /i = 37 
T) will raise the number of as well as \to) and finally 
suppresses the magnetization ruz as shown in Fig. 12b). 

In Fig. |2Ic) we plot the field dependence of niz ob- 
tained by both QMC and MFT, which increases mono- 
tonically as h increases, consistent to experimental re- 
sults. While TJizih) is more or less linear in the mean field 
case, its slope reduces and then increases in the QMC cal- 
culation. This change of slope is again a consequence of 
the formation of condensate. Without the condensation, 
mz{h) will be simply a straight line at T = 0. But the for- 
mation of condensate below Tc increases niz for h < hm 
as discussed above. While the minimum of mz{T) dis- 
appears at he {Tc — 0) and hm, there is no increase of 
magnetization due the condensation at these two fields. 
Therefore the resulted mz{h) curve for h < h„i must be 
convex as shown. Similar argument applies to the high 
field regime h > km and a concave curve is predicted. 
The slope of mz{h), that is the susceptibility, is shown 
in Fig. 131 The small asymmetry indicates the existence 
of t_ states which becomes more significant for smaller 
intra-dimer coupling J as will be discussed in the next 
section. 

Also plot in Fig. \^c) is nt+ and n^, the average num- 
ber of triplet i+ and singlet s per dimer respectively. In 
the MFT, nt+ — v'^p, rit- = v'^g'^ and Hg = u^, where 
M, V, f and g are parameters given by the condensate 
ground state |5'o) (Eq. [S}. rit- is found to be negligi- 
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FIG. 4: (color online) Superfluid density no, which is equiva- 
lent to m^y in the MFT, as (a) a function of temperature and 
(b) a function of magnetic fields. 

ble and therefore is not shown in Fig. ^c). While Ug 
reduces as increases, their total number is almost 1 
(less than 2 percent in difference) for all fields. The in- 
tercept of both curves gives an estimate of hm — 38.6 T 
in which nt+ w rig « 0.5. 

The fact that nt+ essentially coincides with means 
the average number of t- is almost zero {g^ is very small) 
in BaCuSi206. This explains why the hard-core boson 
model b gives results consistent to experimental data. 
However, the derivation from the hard-core boson model 
will be more apparent when we discuss the ruxy where 
the difference is proportional to g instead of g^. 

Temperature dependence of ng is shown in Fig. ^ ng, 
equivalent to the superfluid density m^y in MF, starts 
to rise at about the same Tc obtained from specific heat 
and spin susceptibility and therefore justifies the idea of 
BEC of for the transition observed in BaCuSi206. The 
broken rotational symmetry below Tc that generates rUxy 
is identical to the broken gauge invariance of \^o)- Also 
shown in Fig. ^ is the field dependence of m^y at low 
temperature. The data shows that although MF calcula- 
tion overestimates the order parameter near the field hm, 
in general it agrees well with QMC result. The predicted 
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FIG. 5: (color online) Results of reduced inter-plane coupling 
J = 2.436 meV. (a) Temperature dependence of and num- 
ber density of s, t+ and i_ states, (b) Phase diagram of no{h). 



MF values of critical and saturation fields are given by: 



AiB.9o^c = VJ(J - 4Ji - 2J2) 
M-B5o^s = J + 4Ji-f2J2. 

Since quantum fluctuation is totally surpassed in the fully 
polarized phase, mean field theory predicts the exact sat- 
urated field hg as expected. Below hg, quantum fluctu- 
ation stabilizes the disordered state and so reduces the 
order parameter away from the MF value. The small dis- 
crepancy around he will become more significant when 
we lower the inter-plane coupling J where MF gives a 
much smaller he than the true value. It is also obvious 
that m'^y is asymmetric about hm, contrast to what is 
predicted by the hard-core boson model that retains a 
particle-hole symmetry. This indicates that t_ states, 
although the average number is small, does play a mea- 
surable role in the ordered phase. 



that significant mixing of is expected and new fea- 
tures are possible. Therefore we extend our study to 
smaller intra-dimer coupling J on the same bilayer lattice 
as in BaCuSi206. One possibility to reduce J experimen- 
tally is the substitution of Si by atoms of larger radius 
but the same chemical valence, in order to enlarge the 
intra- dimer distance. Germanium, the one right below 
Si in the period table, is the most natural candidate. 

Reducing J also implies the reduction of he because 
triplet states have lower energies with small J. We choose 
the reduced J to be J = 2.436 meV such that the critical 
field he predicted by MFT is exactly zero. The calculated 
curve (Fig. EJ) of m,{h) of MFT and QMC lay close to 
other each although QMC yields a finite he instead of 
zero. The average number of i_ is not negligible any more 
in this case, no, or rrixy^ shows a much larger discrepancy 
between QMC and MFT than the large coupling J case. 
MFT overestimates rrixy for a large range of magnetic 
field and only when close to the saturation field hg that 
MFT is close to exact results. The asymmetry of rrixy is 
also enhanced when J is reduced as expected. Although 
?i_ is still very small that ruz deviates only slightly from 
n+, the square root of it, ^/nZ = vg, can be as large as 
0.14 and hence implies that Hi, is not valid any more. 
The change of slope of mz{T) is also observable but is 
significantly asymmetric around hm ~ 20 T as shown in 
Fig. El 

Here we discuss how our results of reduced J can 
be successfully applied to a related compound NiCl2- 
4SC(NH2)2 (DTN)^-^ Although this material is not a 
spin dimer system, the single ion anisotropy {D ^ 0.88 
meV) plays the role J and splits the Ni 5 = 1 spin state 
into Sz = Q and = ±1 states, which corresponds to 
s and t± states in the above discussion. Under an ex- 
ternal magnetic field the same field induced ordering has 
recently been observed experimentally in DTN^^. In this 
case, our proposed condensate wave function becomes ex- 
act due to the absence of a corresponding state. There- 
fore the field induced transition observed in DTN shares 
the same physics as in spin dimer systems just described. 
Since D ^ 0.88 meV is comparable to the average ex- 
change coupling J' 0.66 meV, a large asymmetry in 
the phase diagram is expected as we considered in the 
dimer system. Recent experimental data indeed showed 
strong asymmetry in Te{h) as well as in the change of 
slope of magnetization curve. Instead of the proposed lat- 
tice expansion in the compound^, the anomalous change 
of slope in magnetization is in fact a consequence of the 
formation of superfluid condensate which either enhances 
or reduces according to the applied field. 



VI. REDUCED EXCHANGE COUPLING J 



VII. SUMMARY 



In BaCuSi206j since Ji/J and J2/J are relatively 
small, the mixing of state is not significant, as least 
for h ^ he- But there are other dimer systems, like 
TlCuCls, where inter-dimer coupling is strong enough 



Working on the original antiferromagnetic Heisenberg 
model, our QMC calculations reproduce the experimen- 
tal finding of the field induced ordering observed in both 
specific heat and spin susceptibility. The calculated Cy 
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shows a A-like transition at Tc with the ampUtude maxi- 
mizes at /i ~ hm and reduces to zero as h approaches he 
or hs ■ This phenomenon can also be described as a Bose- 
Einstein condensation of magnons. Instead of hard-core 
boson, we propose a semi-hard-core boson, for which up 
to two occupancy is allowed, is responsible for the field 
induced condensation. We showed that the superfluid 
density no is identical to the in-plane staggered magne- 
tization ra^y, which is the order parameter in the spin 
language. Our QMC results show a rise of no = m^y at 
the same Tc of Cv and justify the idea of Bose-Einstein 
condensation. Due to the smallness of Ji and J2 in 
BaSiCu206, the number of i_ states is tiny and leads to 
the success of hard-core boson model Hb. Nevertheless, 
the inadequacy of Hi, is still observable from the asym- 
metry of the phase diagram. A bond operator mean field 
approach agrees well with QMC data of rriz for all field, 
and of for magnetic field close to he and hg. Quan- 
tum fluctuations that ignored in the MET lead to an 
overestimated m^j, around hm however. The difference 
is more significant as one considers a smaller inter-plane 
coupling J. MET predicts a much smaller he than the 
value found by the QMC. The asymmetry of rn^y is also 
much enhanced and signals the failure of simple hard-core 
boson model. This is indeed the case observed in NiCl2- 
4SC(NH2)2- The condensation of magnons also leads to 
the change of slope in mz{h) which is again symmetric 



around hm when J is large but becomes asymmetric for 
smaller J systems such as observed in NiCl2-4SC(NH2)2- 
In the ordered phase, the bond operator mean field 
theory successfully account for longitudinal fluctuation, 
which is absent in the conventional mean field theory, 
in addition to the transverse fluctuation. The longitudi- 
nal fluctuation, or the spin amplitude mode, mixes with 
the transverse fluctuation, the phase mode, at lower field 
but becomes dispersionless and separated from the phase 
mode at higher field. However, topological excitation, 
corresponds to vortices of the condensate, can not be de- 
scribed in the framework of bond operator MET and fur- 
ther work will be needed. Recently a BEC of magnons in 
frustrated triangular lattice of antiferromagnet has been 
proposed and a similar A-transition in specific heat is also 
observedi^. It raises the question whether the BEC is a 
more general phenomenon of quantum antiferromagnet in 
external field and also raises the possibility of some novel 
phenomena arising from quantum nature of the conden- 
sate. 
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